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Abstract 

Dark matter kinetic decoupling involves elastic scattering of dark matter off of leptons and quarks 



^ ' in the early universe, the same process relevant for direct detection and for the capture rate of dark 

matter in celestial bodies; the resulting size of the smallest dark matter collapsed structures should 
thus correlate with quantities connected with direct detection rates and with the flux of high-energy 

r~| ' neutrinos from dark matter annihilation in the Sun or in the Earth. In this paper we address 

' ' this general question in the context of two widely studied and paradigmatic weakly-interacting 

P-c 

(— I , particle dark matter models: the lightest neutralino of the minimal supersymmetric extension of 

the Standard Model, and the lightest Kaluza-Klein particle of Universal Extra Dimensions (UED). 

^ • We argue and show that while the scalar neutralino-nucleon cross section correlates poorly with 

^~^ ■ the kinetic decoupling temperature, the spin-dependent cross section exhibits a strong correlation 

in a wide range of models. In UED models the correlation is present for both cross sections, and 

^^ ' is extraordinarily tight for the spin-dependent case. A strong correlation is also found, for both 



models, for the flux of neutrinos from the Sun, especially for fluxes large enough to be at potentially 
detectable levels. We provide analytic guidance and formulae that illustrate our findings. 
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I. INTRODUCTION 

Weakly interacting massive particles, or WIMPs, have become paradigmatic in the con- 
struction of models for the particle nature of dark matter. Particles with mass in the 
hundreds of GeV to few TeV range, and interacting via Standard Model weak interactions, 
can naturally have a thermal relic density in a range that includes the observed cosmological 
density of the mysterious dark matter. This derives from a thermal history where WIMPs 
were once in thermal equilibrium with the high-density and high-temperature primordial 
plasma - a condition dependent upon the pair-annihilation rate F ~ ^x*^^ being much 
larger than the Hubble expansion rate H] as the temperature dropped below a fraction of 
the mass of the WIMP, the equilibrium number density decayed exponentially with temper- 
ature, as dictated by its Maxwell-Boltzmann equilibrium distribution. Shortly thereafter, 
the precipitous decline of F brought it below if, causing the freeze-out of the WIMPs. The 
resulting number density today is then a function of a combination of effective couplings 
and masses such that, for WIMPs^, one obtains a relic abundance parametrically close to 
the observed dark matter density. 

What described in the previous paragraph pertains to the so-called c/iem«ca/ decoupling of 
WIMPs: the WIMP number density n^ ceases to follow the equilibrium distribution once the 
pair-annihilation and pair-creation rates go out of equilibrium (i.e. they occur less frequently 
than once per Hubble time). After chemical decoupling, WIMPs do not entirely forget about 
the surrounding thermal environment: elastic scattering processes where a WIMP scatters 
off of, for example, a light lepton / (%/ ^ %/) keep WIMPs in kinetic equilibrium. WIMPs 
continue to trace the thermal background kinetically, and structures cannot start to form 
via gravitational collapse. When the rate for elastic scattering processes also falls out of 
equilibrium, structures eventually start forming, and a small-scale cutoff is imprinted in the 
power-spectrum of density fluctuations in the universe. This cutoff scale also defines the 
size of the smallest possible dark matter halos ( "protohalos" ) , some of which might survive 
and populate the late universe, with potentially important irnplications J2]- 

Kinetic decoupling of WIMPs was first discussed in Ref. |3| for heavy neutrinos as dark 
matter candidates, and for supersymmetric neutralinos, first in Ref. j^ some time later. It 



^ but, evidently, also for other particle models with a similar combination of masses and couplings, generi- 
cally also known as WIMP less dark matter models [l|. 



was subsequently argued in Ref. (5| that the typical kinetic freeze-out temperature could be 
as low as a keV, a value that would yield a cutoff scale on the same order of the mass of 
dwarf galaxies - the smallest observed dark matter halos. This would have been a profound 
result, potentially impacting our understanding of the mismatch between the predicted 
and observed number of small-scale dark matter halos in cold dark matter cosmology 6|. 
Unfortunately, Ref. [2] pointed out important kinematic effects that were neglected in 5|], 
leading the latter analysis to vastly overestimating the cross sections relevant for kinetic 
decoupling. The kinetic decoupling temperature calculated in [7| pointed, instead, to the 
MeV to GeV range, with a resulting cutoff scale significantly smaller than dwarf galaxies 
halos, and on the order of the Sun's mass or small fractions of it. 

A number of more recent studies addressed the question of calculating the kinetic decou- 
pling temperature with increasingly finer detail, see e.g. Ref. [8|, including WIMP models 
beyond supersymmetric neutralinos J9| as well as addressing the question of how to connect 
the kinetic decoupling temperature to the scale at which the matter power spectrum is effec- 
tively cut off |10l-ll3|. These studies were paralleled by a series of N-body simulations that 



targeted the nature and fate of the smallest dark matter halos, starting with Ref. 1J| and 



continuing in 



15 



17] . Further analyses studied the question of whether the smallest- scale 



halos would survive tidal stripping and stellar encounters, and whether they would then be 
potentially hovering around in today's galaxies 17l |. with potentially important implications 



for indirect 18l] as well as for direct [19| dark matter detection. Other work also targeted 
the direct detection of these primordial halos (alternately named protohalos, mini-halos or 
micro-halos: the latter two names allude to the size of the halos, which str ong ly depends on 
the particle physics model, see e.g. J9|) via gravitational lensing (e.g. 20l. l2l|). 

In the present analysis we point out that there might be orthogonal handles to pinpoint 
the size of primordial dark matter halos, and thus of the effective cutoff scale of structure 
in the universe. Our main observation is that the same class of processes entering kinetic 
decoupling - namely, elastic scattering off of light fermions - also enters the cross section 
for direct dark matter detection, which is determined by elastic scattering off of quarks 
inside nucleons. Additionally, in a situation of capture-annihilation equilibrium, the rate 
of high-energy neutrinos expected from the capture and annihilation of WIMPs in celestial 
bodies such as the Earth or the Sun also depends on the same scattering cross section. We 
therefore ask, in the present study, whether the mass of dark matter protohalos correlates 



with quantities that could be measured by experiments on Earth, be it via direct detection 
or with neutrino telescopes. 

Here, we take a model-dependent view, and focus on two specific and well-defined WIMP 
scenarios: the lightest neutralino of the minimal supersymmetric extension to the Standard 
Model (MSSM) [22i], an d the lightest Kaluza-Klein (KK) excitation of Universal Extra Di- 



mensions [see Ref. 



58 



% for a review] (We will take a model-independent look at the same 
problem, based on an effective field-theoretic setup, in a forthcoming study [23]). For these 
two paradigmatic WIMP setups, we study correlations between the dark matter cutoff scale 
and rates for direct and indirect dark matter detection. We find that strong correlations 
exist for some quantities, and not for others. We discuss approximations and analytical for- 
mulae that help understanding our detailed numerical results, and we conclude that "earthly 
probes" of the size of the smallest dark matter halos are, in principle and with the model 
assumptions we detail here, possible. 

The ensuing study is articulated as follows: we review in sec. Ill Bl the calculation of both 
the kinetic decoupling temperature and the cutoff scale as a function of this temperature, 
along with potential ways to directly measure the size of the smallest collapsed dark matter 
structures; sec. |TTT] and [IV] discuss in detail, respectively, the case of supersymmetric neu- 
tralinos and of KK dark matter; finally, we discuss our results and draw our conclusions in 
sec.lVl 

II. THE FORMATION AND DETECTION OF PROTOHALOS 

A. Temperature of kinetic decoupling 

The most thorough and comprehensive method to calculate the temperature of kinetic 
decoupling Tkd is a numerical one described in Ref. J2|. This treatment begins with the 
Boltzmann equation in a fiat Friedmann-Robertson- Walker spacetime: 

E{dt-Hp-Vp)f = C[f]. (1) 

Here / is the WIMP phase-space density, E and p are the WIMP energy and comoving 
momenta respectively, and H is the Hubble parameter. C[f] is the collision term, and to 
find Tkd, the necessary C[f] is that of the scattering of a massive WIMP off of a standard 
model (SM) particle that is in thermal equilibrium with the plasma in the early universe. In 



[2], to lowest order in p^/-E^ and the SM particle momentum, this collision term is shown 
to be of the form 

C[f] = c{T)Ml [M^TVl + p • Vp + 3] /(p). (2) 

c{T) is an expression which contains the scattering amplitude of the WIMPs off all possible 
SM scattering partners. We can define a temperature parameter 

^X= o(7^> = T7iT— / 77^P'/(P)- (3) 



3 \2M^/ SM^n^J (27r)3' 
Before kinetic decoupling, WIMPs are in thermal equilibrium with the heat bath, and there- 
fore T^ = T. After kinetic decoupling, the rate of WIMP scattering off SM particles drops 
below the level which is needed to keep them in thermal equilibrium, and so the WIMPs 
cool down due to Hubble expansion, with T^ oc T'^/M^ oc a~^. The transition between these 
two asymptotic behaviors is rapid and corresponds to the temperature of kinetic decoupling, 
Tkd- To find when this change occurs, Eq. [1] is multiplied by p'^/E and integrated over p. 
Using integration by parts, this can be shown to give an equation describing the evolution 
of T^ with the temperature of the universe: 

(9, + 5i7)r^ = 2M^c(T)(T-Tj. (4) 



The author of [2|] has developed a routine which interfaces with the DarkSUSY code [24 1 
and numerically solves this equation. By equating the limiting behavior of T^ in the two 
regimes described above, when T^ = T and when T^ oc T"^ jM^, Tkd is found. It is important 
to note that often before kinetic decoupling the universe passes through the the QCD phase 
transition at T^. ^ 170 MeV, during which the number of relativistic degrees of freedom 
decreases substantially and scattering interactions with quarks are suppressed. To deal with 
this, the code only considers scattering off quarks when T > 4Tc, and after this point all 
scattering is assumed to be with leptons. 

It has also been shown in Ref. 13| that an analytic solution for T^j can be found. For this 
to be done, certain simplifying assumptions need to be made; namely, SM scattering partners 
are relativistic, variations in the universe equation of state are ignored (i.e. Qeff, the number 
of relativistic degrees of freedom, is assumed to be constant), and the scattering amplitude 
is of the form \Ai\ oc (w/M^)" where u is the energy of the SM scattering partner. With 
these assumptions, the solution to Eq. |4] is 



z 



l/{n+2) 



"-^ ^ I z=(a/n+2)(T/M;^:)"+2 



T^ = T 1 -—eMAn-{n + 2)-\z\\ , (5) 



where n is the power of the leading term in \Ai\ and a is a term that contains the leading 
coefficient of the scattering amplitude and Qes- In the limit T — > 0, this equation becomes 



T 



n + 2 



l/(n+2) 



n+ 1 
n + 2 
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(6) 



Tkd occurs when the above limiting behavior matches the high temperature behavior T^ = T. 
Therefore 



Ti 



kd 



Mv 



n + 2 



l/{n+2) 



n + 1 

n + 2 



(7) 



In our work we will use both the numerical code^ and the analytic approximation of Eq. [7] 
to find Tkd- 



B. Protohalo size 

In the period before kinetic decoupling, WIMPs behave as a fiuid coupled to the cosmic 
heat bath via scattering interactions with standard model particles. This coupling leads to 
bulk and shear viscosity in the WIMP fiuid which damps out the primordial structure in 



the fiuid 
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of the form 



he amount of this damping has been shown to be given by a damping term 



111 
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wimp 



(fc, ^kd) 
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wimp 



(^,^i) 



exp 



(8) 



In the equation above, rj is the conformal time and A„imp a function which quantifies the 
amount of fiuctuation in the WIMP density over an isotropic state, so A„imp(A;, r/i) is the 
initial primordial value of the density perturbation function. The characteristic damping 
scale fcd is given by: 



kd ^ U 



m. 



flkd 



H, 



3.76 X 10^ 



kd 



m^ 



Ti 



kd 



(9) 



T],dJ ao ''"" IMpc VlOOGeV/ VSOMeV^ 

After kinetic decoupling, there are no longer interactions between the WIMPs and the 
cosmic heat bath, so damping no longer occurs due to viscosity in the fiuid. However, in this 
epoch free streaming effects are found to significantly damp out density perturbations. In 
this process, the dark matter particles propagate from regions of high density to low density, 
smoothing out inhomogeneities. By considering the coUisionless Boltzmann equation and 



^ We thank T. Bringmann for providing us with his routines. 



using the results of the viscosity calculation as initial conditions, a characteristic scale for 
this damping process can be found, and this comoving scale approaches a constant value 
after matter-radiation equality |ll| : 



kf, 



"eq 



H, 



eq; 



(10) 



ln(4aeq/akd) cto 

where Ogq is the scale factor at matter radiation equality. The damping term for this free 
streaming D^s is of a similar form as D^, and to find the total damping term, the two are 
multiplied together, i.e. D{k) = Dd{k)D{s{k): 



D{k) 



a: 



3 V^fs 



exp 
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(11) 



'wimpl,'^) Vi) 

Comparing kfs and k^, one finds that kfg ^ k^: it is therefore /cfs that determines where 
the exponential cutoff in the mass spectrum is. Therefore, to find the mass of the smallest 
protohalo allowed by these processes, one just calculates the mass of WIMPs contained in a 
sphere of radius vr/Zcfs, i.e. |2|: 
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-Px 



2.9 X lO'^'Mc:, 



1/4 



1 + In (^(^XTkd/SO MeVj /18.56 
mJlOO GeV)'/'(7eff'(Tkd/50 MeV) 



1/2 



(12) 



In later papers 10|, ll2|], it was shown how an additional damping scale is set by acoustic 
oscillations in the cosmic heat bath itself. The calculation of this damping scale is inde- 
pendent of the other two, k^ and kfs, that we previously presented in equations [9] and [in] 
respectively. These oscillations, which are remnants of the inflationary epoch, couple to 
modes of oscillation in the WIMP fluid with k values large enough that they enter the 
horizon before kinetic decoupling. These modes in the WIMP fluid then oscillate with the 
acoustic modes in the heat bath and are damped out, while modes with k values that cor- 
respond to a distance larger than the horizon size at kinetic decoupling do not experience 
such a damping and grow logarithmically. The damping scale for this process is just the 
size of the horizon at kinetic decoupling [k^o ~ TriJ^d); and therefore the cutoff mass for this 
process is the mass of WIMPs enclosed by the horizon at the kinetic decoupling time |2|: 



M^^ 



3 H^ 



3.4 X IQ-^Mq 



T=n. 



rj, 1/4 

50 MeV 



-3 



(13) 



Depending on the parameters of the WIMP model, either Mao or Mfs can be larger, so to 
find a cutoff mass both are calculated and the larger one is used, i.e. Mcut = max[Mfg, Mao]. 



C. Probes of the mass cutoff scale 



A comprehensive review on the detection of sub-solar-mass dark matter halos is given 



in Ref. 26|. Some of these detection methods might provide a more or less direct way to 
infer the value Mcut, even if most of the studies mentioned here do not claim to probe cutoff 
masses as small as Mcut- Dense, nearby protohalos could host enough dark matter pair- 
annihilation to be visible as gamma-ray sources, as envisioned in a number of studies, e.g. 



27H33|. Small-scale subhalos could also contribute to the local cosmic-ray electron-positron 



population, potentially producing a sizable amount to be relevant 
in the abundances of these cosmic rays at 10-100 GeV energies 



34 



or the reported anomalies 



35| . The possibility that 



gamma-ray data would be able to determine the proper motion of protohalos (not necessarily 



on 



but it was shown in 



y the smallest protohalos, however) was first entertained in Ref. 
37( 1 that the diffuse gamma ray background makes this idea unfeasible in practice. Rather 
than aiming to resolve individual substructures, Ref. |3^ considered the anisotropy in the 
diffuse gamma ray emission, arguing that it could be possible to use a statistical analysis 
to measure the substructure mass function. Recently, direct observational constraints from 
the Fermi LAT Collaboration were reported in Ref. (39| in the form of a search, leading to 
a null result, for unassociated gamma-ray sources with spectra that could be conducive to 
particle dark matter annihilation. 

If the Earth were to pass through a dark matter clump, this would lead to an enhanced 
direct detection rate (which scales directly with the local dark matter density), although 



m 



40| it was shown that the presence of substructure in the Milky Way halo is expected, 
on average, to reduce the direct detection rate relative to the rate with a smooth halo and 
no substructure. A long duration direct detection experiment might in principle detect 



variations in the rate due to intervening substrcture, as envisioned in 26|. 



It has been noted in 21| that substructure should affect pulsar timing measurements. 



with M„,t having an effect on the amount of the frequency shift. It has also been discussed 
in Ref. [20| that when there is a time-variable compact source that that is multiply imaged 
by strong gravitational lensing, small perturbations in the gravitational potential, such as 
those caused by protohalos, can lead to variations in the images which could be used to make 
statements about the size of the protohalos. Nanolensing from sub-solar-size dark matter 



halos was discussed in Ref. J41[ , together with the possibility of detecting events with much 



8 



shorter durations and smaller amplitudes than the microlensing events due to stars with 
future surveys. Note that none of these studies would directly provide a probe of the size of 
the small-scale cutoff in the matter power spectrum. 

III. NEUTRALINO DARK MATTER 

We first consider correlations between direct detection rates and protohalo size for MSSM 
neutralino dark matter. To calculate the dark matter direct detection rates we use the 
routines in the numerical package DarkSUSY (see Ref. 2J]; for an in-depth description of 
the direct detection calculation see also Ref. |42|), while Tkd is calculated numerically as 



described in section IIIAI We define our MSSM models by 9 parameters given at the weak 
scale: /i, Mi, M2, M3, rriA, tan/3, ?7isq. At and Af,. This is the same parameterization as 



the "MSSM- 7" described in J2^ (to which we refer the Reader for further details), with the 
change that we let Mi, M2 and M3 vary freely, while in the MSSM-7 the two parameters 
Ml and M3 are related to M2 through GUT-scale gaugino mass universality relations. The 
parameters /i. Mi, M2, M3, rriA and rrisq are scanned over logarithmically in the range of 50 
GeV to 5 TeV, with M2 and /i allowed to take positive or negative values, tan /3 is scanned 
logarithmically over the range 2 to 50, while At and Ah are scanned over linearly in the range 
of -5 to 5. 

All of the models we present in are checked against the accelerator and other particle 
physics constraints contained in the most recent version of DarkSUSY, 5.0.5. They are also 
checked to see if they satisfy the 5a bounds on the relic density from the most recent seven 
yea r release of WMAP data, in which Q^h is constrained to the values .0840 < fi^/i < .1400 



43( 1 . The relic density for each model is calculated with coannihilations using the routines 



in DarkSUSY [2j]. Current (solid line) and future (dashed line) sensitivities from direct 
detection experiments are also included on many of the the plots. For plots with spin 
independent scattering cross sections, we present the current sensitivity of the XenonlOO 



experiment from 4J] and the projected sensitivity of the XenonlT experiment found at [45 



For spin dependent plots, we present the current sensitivity of the 4kg COUPP detector 



47|. 



and the expected sensitivity of the future 60 kg COUPP detector 

The elastic spin-independent neutralino-nucleon cross section depends, at the microscopic 
(quark) level, and at tree-level in perturbation theory, on two sets of diagrams: (i) Higgs 
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FIG. 1: A scatter plot showing the correlation between the neutralino- proton spin- independent (left 
panels) and spin- dependent (right panels) cross sections with the kinetic decoupling temperature 
(upper panels) and with the cut-off scale mass (lower panels), for a large sample of supersymmetric 
models. See the text for details on the definition of the quantities plotted and for details of the scan 
over the supersymmetric parameter space. All sensitivities presented in this plot are for a 100 GeV 
WIMP. 



exchange 
exchange 



including, in absence of CP violation, the two CP-even Higgses) and (ii) squark 



48 



|. Elastic spin-dependent (axial) interactions are also mediated by squark ex- 
change, as well as by Z exchange. Processes relevant to elastic scattering of neutralinos off 
of light leptons and quarks depend on all scattering processes. Since kinetic decoupling typ- 
ically occurs at low temperatures, where heavier fermions no longer participate in the ther- 
mal bath, Yukawa-suppressed Higgs-exchange processes are generically subdominant with 
respect to Z and squark/slepton exchange. This consideration leads us to anticipate that 
the correlation between kinetic decoupling temperature T^d and the spin-independent elas- 
tic neutralino-proton cross section asi be weaker than the correlation with spin-dependent 
processes, ctsd- 
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The theoretical anticipation is in fact confirmed by the results of the extensive scan 
over MSSM parameters we carried out, shown in Fig. [H The upper panels correlate the 
scalar and axial cross sections with T^d? while the lower panels with Mcnt- As expected, 
although a general trend is present, we do not find a tight correlation for scalar interactions 
(panels to the left), while a correlation is definitely present for axial interactions (panels to 
the right, especially for large and potentially experimentally interesting values of the cross 
section). The scatter in the correlation between Tkd and usd is within a factor 2 down to 
o"SD ~ 10~^° cm^, and grows for smaller values of the cross section. The correlation has 
a very small spread at values of the cross section currently probed by the most sensitive 
detectors (for the most recent results from COUPP, probing cross sections as small as few 
xlO^^^ cm^ see [4(|). We investigate and discuss the origin of this scatter in the following 
subsections. Note that o"sd ~ 10~^° cm^ (for a 100 GeV WIMP) corresponds approximately 
to the projected reach of a large (1 cubic meter) DMTPC detector with 50 keV threshold 



operating for one year 491]. This corresponds to a jump of about three orders of magnitude 



over the current detector performance (this does not include indirect limits from neutrino 
telescopes) [46 



46i,y. 



A. The role of the neutralino and of the squark mass scale 

There are two main sources for the scatter in the correlation found between T^j and asB' 
the neutralino mass and the squark mass scale. ^ We discuss these effects both analytically 
and numerically in this section. 

In the limit of heavy squark masses, kinetic decoupling and elastic neutralino-nucleon 
axial scattering are only mediated by Z-exchange, and should thus be tightly correlated. 
However, at a fixed value of o"sd, corresponding to a fixed value of the neutralino-Z cou- 
pling, Tkd inherits a dependence on the neutralino mass scale beyond that produced by the 
dependence of o"sd on M^, resulting in a scatter in the values of T^d for a given value of o"sd- 
We illustrate this effect (for both Tkd and Mcut) in the left panels of Fig. [21 For this scan, 
we use the set of parameters as before, with the changes that the pseudoscalar Higgs boson 
mass TJiA is set to 1000 GeV, the trilinear couplings At = Af, = 0, and, most importantly. 



Note that we always assume that squarks and sleptons are degenerate in mass. 
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FIG. 2: Correlation between the neutralino-proton spin- dependent scattering cross section and the 
kinetic decoupling temperature (upper panels) or cut-off scale mass (lower panels). The panels to 
the left assume a common, large value for the squark mass rusq = 10 TeV, hence the dominant 
process for both kinetic decoupling and scattering off of protons is via Z exchange. The color coding 
indicates three ranges for the lightest neutralino mass, while the black line indicates the analytic 
formula of Eq. ((J7| j. In the right panels we illustrate the effect of lower squark/ slepton masses, 
with m-sq = 0.5, 1, 5 and 10 TeV corresponding to the four color codes. See the text for further 
details on the supersymmetric parameter space scan procedure. All of the models presented in this 
plot have values of Vl^h within or less than the WMAP-7 5a range. The experimental sensitivities 
are for a 100 GeV WIMP. 

m-sq is set to the high value of 10 TeV. The color-coded dots show models with neutralino 
masses in the M^ < 100 GeV (red), 100 GeV < M^ < 500 GeV (orange) and M^ > 500 
GeV (black). 

The dependence of T^^ in the limit of heavy squarks can be understood analytically 
when we calculate both the neutralino proton scattering cross section and Tkd using just 
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the Z exchange tree level diagram, ignoring contributions from slepton and squark exchange 
diagrams. The cross section of scattering a neutralino off a proton when Z exchange is the 



only diagram is given by (following, e.g. |48|): 



3 {gziifMl ' 



^XP,SD = -^ ^4 " E (^L - 9zJ\ ■ (14) 

^ \q=u,d,s / 

To analytically approximate the kinetic decoupling temperature, we use the prescription 
described in section III Al which gives us 

— ^((i)"K!))"- '^^) 



where 

a 



31 /Stt^ jgziifM ^Mpi 



X 
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/ 



E^^*^(fe/)' + fe/)T (1^) 



Here / is all possible SM scattering partners and f^sM is the number of degrees of freedom 
for the SM partners. Combining these expressions, we find 

_ 1 /252 H^ Mi £,((gi„-g|J'A.) M^ \ 

The Tkd oc (^x/'^sd)^^'^ behavior can be seen in the left hand side plots in Figure |2l 
To illustrate the validity of this approximation, we plot this expression as the black line 
in Figure 121 with the sum over fermions including all SM leptons except the quarks, and 

1/2 

5'eff = 4, a value which corresponds to T^d ~ 70 MeV. M^ is set to 100 GeV, and from the 
upper left hand plot, we see that the analytic approximation follows the numerical result 
well for low T^d; with the validity of the approximation becoming less valid at high T^d 
because of that fact that before the QCD phase transition there is quark scattering which 
we do not consider and (?eff varies significantly from the value we chose. 

As squark and slepton masses are lowered, the relative contribution from squark, but 
especially slepton exchange in the kinetic decoupling process increases. As a result, for suffi- 
ciently low squark/slepton masses (which we assume to be at the same scale) the correlation 
between T^d and usd is lost. Namely, we expect T^ to be driven by the sfermion mass, 
while o"sD is tuned by the Z-neutralino coupling and is relatively insensitive to the squark 
mass scale. We illustrate this numerically, and we assess the importance of this factor in a 
robust determination of Tkd from a measurement of usd, in the right panels of Fig. [21 The 
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color-coding shows models where we fix the sfermion mass scale to 10 TeV (red), 5 TeV 
(yellow), 1 TeV (blue) and 500 GeV (green). The figure shows that deviations from the 
expected correlation arise at asD ~ 10""^^ for 5 TeV sfermions, and at asD ~ 10^"^^ for 1 
TeV sfermions. For sub-TeV sfermions the correlation can become weaker, although given 



the negative results on squark and slepton searches from the LHC [50(] we expect little effect 
from light sfermion contributions if usd is close to the range that could be probed by next 
generation experiments, ctsd ~ 10^"^°. 

We also consider the situation in which there is no Z-exchange in scattering processes, 
corresponding to the limit of a purely bino-like neutralino. To obtain models with such 
composition, we use a modified version of the original parameter set where we enforce 
H = M2 = 10 TeV, rriyi = 1 TeV, and At = Af, = 0. This leaves us with the parameters Mi, 
M3 and rrisq, which are scanned logarithmically over the range 50 GeV to 5 TeV, and tan/3, 
which is scanned logarithmically over the range 2 to 50. A plot of Tkd and Mcnt versus ctsd 
for these models is shown in Fig. [3l As binolike neutralinos are very weakly interacting, 
only about 100 of the 2.5 x 10^ models in the plot do not produce a relic density greater 
than the WMAP-7 5a range. 

For such a class of models, we see a correlation between T^d, Mcut, and aso, with, just 
as before, an additional dependence on the neutralino mass. Unlike in the high rrisq case, 
there is a scatter of large neutralino mass models down into the smaller mass bands. We 
have identified these models as having a small splitting between M^ and nisq. In [13|, it was 
shown that for binolike neutralinos, T^d is of the form: 

T,, = 7.5MeV -^-1 77^^7.7 (18) 



Ml I \ 100 GeV J 

Since in the MSSM-7 parameterization the squark mass scale is the same as the slepton 
mass scale, when [m'^„ — M^) /M^ <^ 1, T^d is driven to a much smaller value than when 
the splitting is large, as can be appreciated from Eq. ( 1T8|) . When Z-exchange is the only 
process relevant for the interaction, the mass splitting does not enter into T^d; and we see 
no similar effect. 

For these models, the highest possible spin dependent cross section is significantly smaller 
than that of the models from the full parameter space scan as well as that for models where 
sfermion exchange diagrams are suppressed. As such, these are not models that the current 
generation of direct detection experiments would explore: this additional source of scatter 
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FIG. 3: A scatter plot showing the correlation between the neutralino-proton spin- dependent scat- 
tering cross section and the kinetic decoupling temperature (upper panel) or cut-off scale mass 
(lower panel) for super symmetric models with bino-like lightest neutralinos, hence with suppressed 
coupling to the Z. The three colors correspond to different ranges of neutralino mass. 

in the correlation between protohalo size and scattering cross section is thus not worrisome, 
at a practical level. Furthermore, for general sets of models where one has both Z and 
sfermion exchange scattering diagrams, these results show that when there is a relatively 
large scattering cross section, the sfermion exchange contribution to that cross section is 
subdominant to the contribution from Z exchange: for supersymmetric models that might 
be detectable with current or future generation detectors it is thus valid to approximate the 
scattering cross section as being due solely to Z exchange. 

B. Neutrino Telescopes 

Spin-dependent neutralino-nucleon interactions drive, quantitatively more than spin- 
independent interactions, the capture rate of neutralinos in the Sun. If the neutralino 
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FIG. 4: A scatter plot showing the correlation between the flux of neutrinos from the Sun from dark 
matter annihilation with the cut-off scale mass (lower panels), for a large sample of supersymmetric 
models. The panels to the right focus on models with large squark masses. The upper panels use a 
threshold of 10 GeV for the neutrino energy, while the lower ones use 100 GeV. See the text for 
details on the definition of the quantities plotted, and for details of the scan over the supersymmetric 
parameter space. 



pair-annihilation rate is large enough so that capture and pair-annihilation in the Sun are in 
equilibrium (which is typically the case across the MSSM parameter space 5l|), the actual 
rate of neutralino annihilation is governed by the capture rate. We thus expect a correlation 
of the rate of high-energy neutrinos resulting from neutralino annihilation inside the Sun 
and Tkd. Note that effects such as the detector energy threshold are expected to impact the 
correlation and to potentially disrupt it. 

We investigate in Fig. H] the correlation between the flux of neutrinos from the Sun 
integrated above two representative energy threshold, namely 10 GeV (upper panels) and 
100 GeV (lower panels). Models in the two panels on the right assume heavy squark masses, 
while those on the left scan over the general MSSM parameter space as before. To calculate 
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the neutrino flux, the routines from DarkSUSY described in 23, l52|] are used. We note that 
a rather tight correlation exists between the neutrino flux from the Sun, (p^, integrated above 
10 GeV and M^ut, as long as (p^ > 10^^ km^^yr"^. The most recent IceCube results looking 
for neutrinos from solar WIMP annihilation through the W^W~ channel claim a sensitivity 
;o a total muon flux from this annihilation of about 3 x 10^ km^ yr^^ for a 1 TeV WIMP 
53| . Using the routines from DarkSUSY, we find that this muon flux corresponds to a range 



of incoming total neutrino fluxes from about 10^° — 10^^ km~ yr~^ for MSSM models when 
the neutrino threshold energy is 10 GeV. 

Larger energy thresholds tend to loose the desired correlation: a very significant depen- 
dence on the neutralino mass is present in the flux above 100 GeV, as, for example, almost 
no neutrinos with those energies are produced for neutralinos with masses below 200 GeV 
or so. Therefore, a model with a given Mcut can well have a vanishing neutrino flux if the 
neutralino mass is small enough! Interestingly, with the deployment of the DeepCore detec- 
tor 5^ the effective energy threshold for neutrino detection of the IceCube system has been 
significantly lowered. Again, for fluxes large enough to be above potentially detectable lev- 
els, we find that a tight correlation with M^ut is present. The correlation we find is expected 
to further improve should plans to deploy an additional, even more thickly instrumented 
section of the detector, PINGU, come to fruition |55|]. 



IV. UNIVERSAL EXTRA DIMENSIONS 



The Universal Extra Dimensions (UEDHramework offers an interesting setup for a WIMP 



model alternative to supersymmetry 56 



57l | (see also Ref. [58[ for a review). The lightest 



Kaluza-Klein (KK) n = 1 excitation, usually the first KK mode of the hypercharge gauge 
boson B^^\ is stable by virtue of the so-called KK parity 58|. The lightest KK particle, 
or LKP, makes for a phenomenologically viable WIMP dark matter candidate. The parti- 



cle properties of the 
will consider here 



iKP^ depend, for the minimal version of the UED scenario that we 



59| . upon three parameters: the effective cut-off scale A, the inverse 
compactification radius 1/-R, and the value of the Standard Model Higgs mass. The latter 
quantity is especially crucial for the calculation of the spin-independent LKP-nucleon scat- 



We shall use LKP and B^^-' interchangeably in what follows. 
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FIG. 5: A scatter plot showing the correlation between the B^^' -proton spin- dependent cross sec- 
tion with the kinetic decoupling temperature (upper panel) and with the cut-off scale mass (lower 
panel), for a large sample of universal extra- dimensional models, and the analytic approximation 
of Eq. [21\) . Plotted sensitivities are for a 1 TeV WIMP. See the text for details on the scan over 
the UED parameter space. 

tering cross section. 1/-R, instead, sets the mass scale of the KK levels, including the mass of 
the LKP, while A feeds in the details of the particle spectrum. Here, we consider the range 
500 GeV < l/R < 1400 GeV, 10 < KR < 40 for our scans. We also scan over the values of 



the Higgs mass allowed by current collider constraints |60l-l62| , with a maximum Higgs mass 
of 600 GeV. We then find the relic density for all models using the results of 631] and check 
these against WMAP relic density constraints. 

As for neutralinos, the same diagrams contributing to elastic 5'^^^-nucleon scattering 
contribute to the process of kinetic decoupling. In particular, spin-dependent scattering 
depends upon KK-quark exchange, while spin-independent scattering is primarily driven 
by processes mediated by Higgs exchange. As a result, the general expectation for UED 
is not dissimilar to what we awaited in the context of supersymmetry: a tight correlation 



18 



between spin-dependent elastic processes and kinetic decoupling, and a looser correlation 
for the scalar cross section. Figures and E] accurately confirm these expectation. In Fig. \5\ 
we show, for 5000 minimal UED models, the correlation between the elastic spin-dependent 
cross section and the kinetic decoupling temperature (top panel) and cut-off mass scale 
(lower panel). All of these models have the same Higgs boson mass, 125 GeV, which does 
not effect the calculation of T^d or Mcut, but does have a large effect on the relic density. Here 
Tkd is once again calculated using the numerical method described in section IIIAl while the 



spin dependent cross section is calculated using the approximate formula from Ref . [571, l58| : 



.42 j,lTeV\V0.1 



In the formula above. A, is the mass splitting between the right handed quarks and the 
LKP, A„ = (M (1) — M^(i))/M^{i), with all of the KK quarks taken to have the same mass 
for simplicity. An analytic approximation for T^d was already found in [l3|, which is 

Tkd-3xlO^MeVAy2(^y^', (20) 

with Ae = (Mg{i) — Mq[i))/Mq[i). Putting together equations (IT9|) and (l20l) . we get an 
analytic approximation relating Tkd and cr^cijp^sD- 

A \ 1/2 / T x 5/8 /in-42 2\ 5/16 



There is no dependence on dark matter mass as there was in the SUSY case, rather 
the kinetic decoupling temperature just goes like Ugj-, , leading to the strong correlation 
displayed in Figure [51 The validity of this approximation is shown by plotting it as the 
black line in Figured with Ag = .17 and Ag = .1. For the relatively low values of T^d we 
find for UED models, M^ut is always set by the acoustic oscillation cutoff of equation (1131) . 
In plotting the analytic approximation for M^ut, we use Qcq = 3.5, which corresponds to 
Tkd ~ 30 MeV. 

In the spin-independent case, to find the cross section we use the approximation from 



m 



.o in 46 sATeVV /lOOGeVV ^ ^^ /^l TeV V /0.1\ 
(^bWv si ~ 1-2 X 10"^^ cm=^ — + 0.09 -^ 



(22) 



The results of using this formula are shown in Fig. [6l where were are now scanning over a 
full range of Higgs masses. Here the scattering cross section at a given value of 1/R depends 
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FIG. 6: Left: a scatter plot correlating the LKP mass with the LKP-proton spin-independent 
scattering cross section (upper panel), and correlating the same cross section with the cut-off scale 
mass (lower panel), for a set of models featuring both LHC-allowed and excluded values for the 
Higgs mass. Right: a scatter plot correlating the B^^' -proton spin-independent cross section with 
the kinetic decoupling temperature (upper panel) and with the cut-off scale mass (lower panel) , for a 
large sample of universal extra- dimensional models, and for a range of Higgs masses. In all panels 
except the top left, the plotted sensitivities are for a 1 TeV WIMP. In the right hand panels, all of 
the models produce a relic density that is within or less than the WMAP-7 constraints on Q^h. 

sensitively on the value of the Standard Model Higgs mass M^, as shown in the top panel to 
the left. In the figure, the green points correspond to values allowed before the most recent 
LHC results on searches for the Standard Model Higgs [60], which are now ruled out (the 
region without points corresponds to values of M^ already ruled out by searches with the 



Tevatron and LEP 



61 



62|). LHC results have therefore severely constrained the prediction 



for the scalar S*^^^-nucleon scattering cross section to a range of about an order of magnitude 
around 10~^® cm^, along with significantly tightening the correlation between 1/R and this 
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cross section. This, in turn, implies a correlation of the latter cross section with the cutoff 
scale, illustrated in the lower-left panel. 

In the future, more accurate measurements of M^ will yield an increasingly tighter cor- 
relation between the spin-independent cross section and the kinetic decoupling temperature 
and cutoff scales, as we show in the right panels. There, we show with different colors models 
corresponding to the ranges 114.4 < Mh/GeV < 120 (cyan), 120 < Mh/GeV < 130 (brown) 
and 130 < Mh/GeV < 141 (red). We also calculate the expected analytic form for the 
correlation between the quantities shown in the right panels. As Higgs exchange dominates 
over the KK quark exchange processes in the calculation of the spin independent scattering 
cross section, we approximate Eq. (!22|) as 

1 TeVV /lOOGeVy 



a„..,,, « 1.2 X 10-« cm= ^^^j ^^j— j . (23) 

Combining equations (!20|) and (123!) . we find 



T„.34MeV^^J ^^^ J ^—— J . (24) 

This approximation for T^j and the corresponding result for M^nt are plotted in the right 
hand side of Fig. |6l for Mh = 125 GeV and Ag = .01. The analytic approximation 
underestimates Tkd slightly, which is due to ignoring the KK quark exchange processes in 
the expression we use for asi- However, the behavior of the numerical and analytic results 
is the same, with T^j oc (M^asi)"^/^ and with a scatter occurring due to the varying mass 
splittings. 

Finally, we note that we estimated the neutrino flux from the Sun in the case of LKP 
dark matter, with results that mirror the same behavior and correlation as we found in the 



case of supersymmetric models, illustrated in Fig. HI It is shown in Ref. pm\ that the event 
rate in neutrino telescopes correlates strongly with M^a) . We have shown that for both the 
spin dependent and independent cases in UED there is also a correlation between M^(i) and 
the scattering cross section. This is shown analytically in Eq. f lT9|) for the spin dependent 
case and numerically in the upper left hand panel of Fig. [6] for the spin independent case. 
Therefore, there is a resulting correlation between the neutrino flux and o"sd or crsi. The 
larger mass of LKP with respect to the possibly light neutralino case produces a slightly 
smaller spread as the one shown in Fig. HJ especially with a larger neutrino energy threshold. 
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V. DISCUSSION AND CONCLUSIONS 

We addressed the possibility of establishing the small-scale cutoff of the cosmological 
matter power spectrum in a variety of particle dark matter models via dark matter direct 
and indirect detection experiments. We argued, and showed with analytical calculations 
and numerical results, that the kinetic decoupling temperature, which sets the cutoff scale, 
correlates tightly with the spin-dependent elastic scattering cross section of WIMPs off of 
nucleons. There also is a generically tight correlation with the flux of high-energy neutrinos 
from the Sun - whose intensity depends on the capture rate of WIMPs in the Sun, in turn 
set by the same axial scattering cross section. A weaker correlation is found in the case of 
scalar WIMP-nucleon interaction. Control over the spectrum and properties of the Higgs 
sector will dramatically improve this latter correlation, especially in the case of minimal 
Universal Extra Dimensions. All correlations we found are tightest when the detection rates 
are largest. 

In summary, we find that "earthly" probes of the small-scale cutoff to the cold dark 
matter power spectrum are possible in the foreseeable future. Multiple measurements, for 
example of both scalar and axial scattering cross sections and of a flux of neutrinos from 
the Sun, would help pinpoint a concordance dark matter model and, in view of the results 
reported here, its cosmological bearings on structure formation. 
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